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SUMMARY 

This report documents a computer program called Ensemble Averaging of Acoustic 
Data, The program samples analog dita , analyzes the data, and displays them in 
the time and frequency domains. Ha^d conies of the displays are the program's 
cutout. The documentation includes a description of the program and detailed 
user instructions for the program. This software was develooed *or use on th 
Ames 40x30- 'oot wind tunnel's Dynamic Analysis System consisting of a PDP-11/ 
computer, two RK05 disk drives, a Tektronix 511 keyboard/display terminal, an 
FPE-- Fourier Processing Element, and an analog-to-di gi tal converter. 

This software was developed by Beam Engineering, Inc. under ourchase order no. 
A33Gc55 to NASA-Ames Research Center witn Marianne Mosher as contract monitor. 





LIST OF SYMBOLS 


DUj) 

F(..) 

HU.) 


i 

K 

KMAX 

N 

n 

P(n) 

S(n) 

S^(n) 

SPL 

VCAL 


Array of data in frequency domain 

Array of power spectrum data after bandwidth has been 
corrected 

Array of data from single-sided transform with correct 
ampl i tude 

Complex conjugate of array 

Array of power spectrum data 

Index for frequency 

Sample index 

Number of averages 

Number of samples 

Index for time array 

Array of data in pressure units of N/m^ 

Average sequencer of sampled time data 
One sequence of sampled time data 
Sound pressure level 

Calibration constant which converts volts to N/m^ 
Frequency 
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INTRODUCTION 


The program "Ensemble Averaging of Acoustic Data" samples analog acoustic data, 
analyzes it, and displays it in both the time and frequency domains. Three graphs - 
a time history graph, a power spectrum graph, and a sound pressure level graph - 
are produced. Hard copies of these displays are the program's output. 

This program was written for use on the Ames 40x80- foot wind tunnel's Dynamic 
Analysis System. The output is written to a Tektronix 611 keyboard/display 
terminal. The data is input from a playback tape recorder connected to the Dy- 
namic Analysis System. This program can also use on-line data. For this report 
the descriptions are given as if data is played back through a tape recorder. 

This report includes a description of the program which explains how the program 
analyzes the data and instructions for operating the program. The program is 
written in Fortran and uses Fortran subroutines for data gathering and processing. 
These subroutines are found in Reference 1. The report also contains a general flow- 
chart for the program in Appendix A, a flowchart of the code in Appendix B, a 
listing of the program in Appendix C, and examples of typical running sequences 
in Appendix D. 

There are two versions of this program: Time Frequency 2 (TF2) and Time Frequency 3 

(TF3). They are identical except for the way the bandwidth is corrected for the 
power spectrum graph and they write the calibration data to different files. The 
differences in TF2 and TF3 will be explained fully in the Program Description 
section and the User Instruction section. This report is valid for both versions. 

I would like to express my thanks to Marianne Mosher for her help and guidance in 
the development of this software. 


PROGRAM DESCRIPTION 

The following is an explanation of how the program "Ensemble Averaging of Acoustic 
Data" analyzes analog acoustic data. The flowchart in Appendix A illustrates the 
process. It would be beneficial to the reader to refer to Appendix A while read- 
ing this section. 

The first step is to calibrate each channel of microphone data. This is accomp- 
lished by first sampling the analog data from a channel specified by Thumbwheel A. 
(Thumbwheel A is on the front panel of the Dynamic Analysis System) The analog 
data are digitized by the Dynamic Analysis System's converters and a sample length 
of 2048 data points are taken. The sampling frequency used is 1024 Hz with a low- 
pass cut-off frequency of 500 Hz. Data acquisition is started on receipt of an 
electronic trigger signal. One average is taken of the calibration data. 

Next, the calibration data are converted from the time domain to the frequency 
domain by a single-sided Direct Fourier Transform. The output values are cal- 
culated by the following formula: 


- 3 - 



D(toj ) = 1/N [P(n)exp[-i2iTa)j n/N] x [v8/3 x 1/2 x (1-cos 2Trn/N)] 

where j - 0 to N/2. The Normalized Hanning operation is performed in the frequency 
domain after the direct transform is computed (Reference 1). 

The correct amplitude from the single-sided transform is then computed by multi- 
plying the data by a factor of 2 

G(-^j) = (2) X (D(.^j)) 

Then the transformed spectrum is converted to a power spectrum by multiplying the 
spectrum with its complex conjugate 

H(cuj ) = [G(u)j )] X [5(ojj )] 

Next, the calibration constant is computed. The calibration data points 400 to 
600, corresponding to the frequencies 400 to 600 Hz, are summed and the square 
root is taken of this sum. This value is then multiplied by the calibration factor 
which is input by the user. The calibration factor is determined by the type of 
calibrator and microphone used. The standard calibration factor is 1.0. The 
value, VCAL, is the calibration constant 


VCAL = Cz^= 4 go^(^j ^ X [calibration factor] 

VCAL is used to compute the pressure constant. This sequence continues until 
all the microphones have been cal i bitted. The program is then ready to sample the 
analog acoustic data. 

The analog data are sampled as in the calibration sequence. The digitizing rate, 
sample frequency, gain, and the number of averages to be used are determined by 
user input. Anti-aliasing filters within the Dynamic Analysis System use a cut- 
off frequency that is determined by the sample frequency chosen. Data acquisition 
starts on receipt of an electronic trigger signal for each sample. 

The data are then averaged in the time domain. The sampled waveforms are added 
up point by point and then divided by the number of averages 

S(n) = Sj^(n)/Kmax 

The data are converted to pressure units of N/m‘ by multiplying the data by the 
computed pressure constant 

P(n) = J ^ 31.7/VCAL) x S(n) 

The analyzed time history is then displayed on the Tektronix screen. 

The data are converted from the time domain to the frequency domain, the amplitude 
is corrected, and the transformed spectrum is converted to a power spectrum by 
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the same methods used in the calibration process. At this point the correct band- 
width is computed. The power spectrum data points are summed together by the fol- 
lowing formula; 

If TF2 is used, i H(o.' 2 j)for j = 1 to n/2 

If TF3 is used, F(a-) = j = 1 to n/3 

where n is the number of frequencies in the spectrum or half the number of points 
sampled. Summing the data points changes the bandwidth by a factor of two if TF2 
is used and by a factor of three if TF3 is used. The power spectrum graph is 
then displayed on the Tektronix screen. 

The sound pressure level graph is displayed on a log scale by the following 
formul a : 


SPL = 101ogio[F(ojj)/(. 00002)^] 

This sequence is repeated for all data to be analyzed until the user exits the 
program. 


USERS' INSTRUCTIONS 

The following is an explanation of how to operate the program "Ensemble Averaging 
of Acoustic Data" on the Dynamic Analysis System. Appendix D contains an example 
of a typical operating sequence. It would be beneficial to the reader to refer 
to Appendix D while reading this section. 

The operator must load the disk pack labeled "TF2 and TF3" onto the RK05 drive and 
boot the Dynamic Analysis System. The operator then initializes the program by 
entering the following command: 

•R TF2 <cr> 

Once this command has been entered the program automatically reads the 
calibration data from the file VCAL2.DAT and stores it in memory. If 
the user wishes to use the TF3 version, R TF3 is entered and the cali- 
bration data are then read from the file VCAL3.DAT. The program will 
then begin prompting the operator for input parameters. The operator 
enters the desired values as described below. 

ENTER TEST NUMBER NNN <cr> 

The operator enters the test number of the analog acoustic data to be 
analyzed. The test number can range from 1 to 999. Entering a zero 
or a carriage return will cause the program to terminate. 
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ENTER RUN NUMBER NNN <or> 

The operator enters the run number of the data to be analyzed. The run 
number entered can range from 1 to 999. Entering a zero or a carriage 
return signals the end of a teH to the program which then prompts for 
the next test number. 

ENTER POINT NUMBER NNN <cr> 

The operator enters the point number of the data to be analyzed. The point 
number can range from 1 to 999. Entering a zero or a carriage return 
will cause the program to prompt the next run number. 

WHAT TYPE OF CALIBRATION? 0 OR N (NO <cr>) 

The operator enters which type of calibration data are used, old or new, 
by entering an 0 or N respectively. No carriage return is necessary. 

If the operator enters an 0, the old calibration data are used and the 
operator is then prompted for the voltage code. If an N is entered, 
the following two prompts appear on the screen. 

ENTER NUMBER OF MICROPHONE TO CALIBRATE <cr> 

The operator sets the appropriate channel on Thumbwheel A and then enters 
the microphone number. ***Before entering a carriage return the user 
turns on the tape recorder so the proper calibration signal is present.*** 
The microphone number can a-ange from 1 to 99. 

WHAT IS THE CALIBRATION FACTOR FOR MICROPHONE XX (REAL INPUT) <cr> 

The operator enters the calibration factor for the microphone just 
calibrated. The standard calibration factor is 1.0. The program then 
prompts for microphone, and calibration factor continues until the 
operator enters a zero or a carriage return in response to the prompt 
for the microphone number. Once this occurs it indicates to the program 
that the calibration process is complete. The new calibration data are 
then written to the appropriate file and the next prompt then appears 
on the screen. 

ENTER VOLTAGE CODE <cr> 

The voltage code specifies the maximum analog input signal level that 
is allowed while sampling the analog data. The voltage chosen should be 
higher then the maximum peak voltage of the signal. The voltage code 
can range from 1 to 7 and corresponds to the volts used as specified 
in the following table. 
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VOLTS 


VOLTAGE CODE 


± .125 1 

± .25 2 

± .5 3 

±1.0 4 

±2.0 5 

±4.0 6 

±8.0 7 


Entering a zero or a carriage return will cause the program to prompt 
for the next point number. 

ENTER SAMPLE FREQUENCY CODE <cr> 

The sample frequency code entered specifies the sampling frequency of 
the analog to digital converter. The sample frequency code can range 
from 1 to 7. The sample frequency code also specifies the cut-off 
frequency to be used by the anti-alising filters according to the follow- 
i ng table. 


SAMPLE 

CUT-OFF 

SAMPLE 

FREQUENCY 

FREQUENCY 

FREQUENCY CODE 

51.2 

20 Hz 

1 

204.8 

100 Hz 

2 

512.0 

200 Hz 

3 

2048.0 

1- KHz 

4 

5120.0 

2 KHz 

5 

20480.0 

10 KHz 

6 

51200.0 

20 KHz 

7 


Entering a zero or a carriage return will cause the program to prompt 
for the voltage code. 

ENTER SAMPLE CODE <cr> 

The sample code entered specifies the number of analog data points to 
be acquired and stored per frame and per channel. The sample code can 
range from 1 to 4 and corresponds to the number of samples taken as 
specified in the following table. 

NUMBER OF SAMPLES SAMPLE CODE 


128 1 

512 2 

1024 3 

2048 4 


Entering a zero or a carriage return will cause the program to prompt 
for the voltage code. 
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ENTER MICROPHONE NN <cr> 


The microphone number entered is the microphone number of the data to 
be graphed. The microphone number can range from 1 to 99. Entering 
a zero or a carriage return will cause the operator to be prompted for 
the number of samples. 

ENTER GAIN (REAL INPUT) <cr> 

The gain must be entered as a real number. Entering a 99.0 will cause 
the program to prompt for the microphone number. 

GAIN IS XXX Y OR N (NO <cr>) 


The operator checks if the correct gain was entered. If the gain entered 
was not correct, the operator enters an N. The program then prompts for 
the gain. If the correct gain was entered, the operator enters a Y 
and the program continues with the following prompt. 

DATA ACQUISITION STARTS WHEN <cr> IS ENTERED 
ENTER NUMBER OF AVERAGES NNN <cr> 

The number of averages entered is the number of averages to be taken 
of the waveforms. ***Before the carriage return is given the operator 
turns on the tape recorder so the proper signal is present.*** The 
number of averages can ran^e from 1 to 999. Entering a zero or a 
carriage return will cause the program to prompt for the gain. 

At this point the time history and spectra graphs are displayed on the Tektronix 
screen. The display will stay on the screen until a carriage return is entered. 
Before a carriage return is entered a hard copy of the display must be made since 
there is no way to back up to the previous display once the carriage return has 
been entered. After the sound pressure level graph is displayed and the operator 
has entered a carriage return the program prompts for the next microphone number. 
The program will continue until a zero or a carriage return is entered in response 
to the prompt for the test number. 

Various error messages may appear on the Tektronix screen while running this 
program. If this occurs, the user should refer to Reference 1. 


REFERENCE 


1. Reference Manual TSALF Time Series Analysis Library - Fortran (for RT-11 
Fortran System), GenRad Time/Data Division, 1976. 
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TYPE 60 

FORMAT ‘ENTER RUN NUMBER NNN 

ACCEPT 40. IRUN 

IF < IRUM.£Q.000> GO TO 20 
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hCCEPT 40.. I Pul NT 
IF < IPuIHT . Eu . 000 > uO TO bti 

CT HEH DR OL.D ChLIBRUTIuH T 
t TTTTTTTT 

TYPE 90 

90 FGRMhKIK. ‘WHhT TYPE OF ChLIBRhTIuH ? 0 OR H (MO <CR>> • 

CALL I PuKE< " 44 , “ I OOOO :• 

100 I CAL = ITTlNRv > 

IF( I CAL . LO . 78 . OR . I CAL . EQ . . '9 > CALL ECHO ( I CAL > 

I F< I CAL . EQ . 78 . OR . I CAL . EQ . 79 > CALL I POKE( " 44 .. “ OOOOO > 

IF( ICAL EQ. 78.' lO TO llO , 

IF( ICAL EQ. 7S‘.' GO TO 140 '' 

GO TO 100 

cm PROMPT OPERATOR FOR HUMBER OF HEXT MICROPHOHE TO CALIBRATE TTY 
110 TYPE 120 

120 FORMAT*: •+ • EHTER HUMBER OF MICROPHOHE TO roLIBRATE ’.T* 
ACCEPT 130.. IMIK; 

130 FORMAT-: 12 ;> 

IF (IMIK.HE.O:* lO to 4Gi0 
ICAL = 79 

C»: HHEH DOME HITH CAL I BRAT I OH HRITE CAL I BRAT I OH DATA TO FILE UCAL2 DAT T 

CALL ASSIGH-: 2.. ‘ UCAL2 . DAT ‘ / 9.. ' HEH * 

HRITE-: 2.. 1 > -: UCAL( I .>.• 1 = 1. 99 > 

CALL CLOSE-: 2 > 



A 


C M E H U \ 

CALL ERASE-: 'KG' > 

01:1:4:4: USER PRONF-TS FOR ‘.'OL.TFiGE CODE FiUD SAMF'LE F'REQUEHCV m 


I 

CM 



1 t‘U rvph 170 

170 FORMAT ’EUTER UUL I hGE CODE ' T • 
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ACCEPT ISO.INDEY'S 

IF (.(.I NDEKS . GE . 1 > . FIND , < I NDE! CS . LE . 7 > > GO TO 2S0 
IF < INDEXS.EQ.O;:- GO TO 160 



I 

o\ 


TYPE 240 

240 FuRMhT < > .. 1S>=: . •NHOrr ShMPLE frequency is 1-7! • > 
GO TO 220 

250 CALL ERASE*: * KB ‘ > 

CKTTT USER PROMPT FOR NUMBER OF SAMPLES TTTT 


260 

270 

1 

I 

2y0 

2yO 


300 


TYPE 270 

FORMAT*: ^ - ..T20, 'NUMBER OF SAMPLES 
. T43,. • 1 • T28.. • 512' . T43 . ' 2 ' /VY 


ENTERS' 
7.. * 1024 • , T43.. 


'2fi4«' ,T43_. '4* > 

TYPE 2y0 

FORMAT < V . lyv, -EINTER SAMPLE CODE ' ,4 > 

ACCEPT 180.. INDEXZ 

IF *: *: inde:^:z . ge i > . and . •: inoekz . le . 4 > > go to 310 

IF •: INDE>;Z .EQ.O> GO TO 140 . 

TYPE 300 

FORMAT •: 18>:.. •MMAT? SAMPLE NUMBER IS 1-41 *;• 

lO to 280 


^ / I 

* ir t 


12 


cm I BANDMIDTH.. MOLTS, AND SAMPLE FREQUENCY OUTPUT TTY 

3 1 0 BDU= SAM*: I NDE7.S I S I ZE( 1 1 !DEXZ > 

TYPE 320. F:DM 

320 FORMAT ( 1 S-K.. ’ BANDH I DTN = * F 1 0 . 3 > 

TYPE 32 1 > UOL.T*: I MDE: 1. > 

321 FORMAT •: - Y ISZ, ' MOLTS = +/- SF5.2;* 

TYPE 322 ,. SAM*: I UDE2S > 

322 FORMAT < - Y ISK . * SAMPLE FREQUENCY = • ..F8.2> 

CYYYYY USER PROMPT FOR MICRuPHuNE, GAIN, AND AMERAGES YYYYYYY 

323 TYPE 324 

324 FORMAT*: -Y lyX,' ENTER NICROPNONE NN 
ACCEPT 180,1 NIK 




‘.J cc» 



A 


r') 

I 

--I 


326 

32? 

329 

329 

33m 


331 


-3-7-7 


ly 


IF (IMIK.EQ.ei:> GO TO 250 
TYPE 326 


FORMhT <■+ ENTER GAIN (REAL INPUT) '.$> 

ACCEPT 327 ..GAIN 
FORMAT (F5.1> 

IF (LAlN.Eu.yA. > LU TO 323 
TYPE 329 .GAIN 


FORM AT( ly.K. 'GAIN IS ‘..F5.1.. ' Y OR N (NO 
C AL.L I POKE( “ 4 4 , " 1 OOOO > 

I DU = ITTINR( 

IF ( I DU . EQ . 7 . 3 . OR . IDU . Eli . 99 > CALL ECHO ( I DU 

IF (I DU . EQ 73 . OR . I DU . EQ . 99 > CALL I POK.E( " 44 

IF < IDiJ.EQ.79> lO TO 325 
IF (IDU.EQ.99> GO TO 331 
GO TO 330 
TYPE 332 

FORMAT-: •-»- * 19>L ‘DATA ACQIJISTION STARTS 

19X. ‘ENTER NUMBER OF AUERAGES NNN ' 

ACCEPT 333.N0AUE 
FORMAT-: 13 > 

IF •' NOhUE EQ . 0 ;• LU TO 325 
lO to 410 


<CR>> •,$> 


/ " 00000 > 


WHEN <CR> 


IS ENTERED'. 


CT DAT A T 

mm CALIBRATION PARAMETERS FOR ACQUISITION OF ANALOG DATA TTTT4; 

400 IHDE2S=9 

INDEX2=4 
INDEXL=5 



Ctm PrtRHMETERS hHD DhTh FOR hCQUI SITIOM OF hNhLOG DhTh 


4 1 0 COL-L hN I N I J< I hF'RhV .. h I E: I , DUfl , DUN .. I S I ZE< I NDE>^Z hi. 

I SFRECK i: r > .. 1 RhHGEv I NDEKS > .. I LEUEL< I NGE>=:L > AER.id.6. 
1 IFILTE-: i:uDE>(S>.. 1 ;> 

I F vIER .Eu .O GO TO 43£i 
TYPE 420.. lER 

420 FDRNhT V . .. 5/: .. ‘ ERROR hW I! N I T ' 1 5 > 

STOP 

430 IF ( IChL.EQ.7S> GO TO 440 

AC COUPLING.. DC OFFSET, AND CALIBRATION OPTIONS FOR 
cm INPUT UF ANALOo DATA 

CALL AN IN IT (2^IER> 

IF ■' IER.EO.O.> LiO TO 440 

TYPE 420.. lER 

STOP 

440 CALL MUXCH < IER> 

IF (I ER NE . 0 > GO T 0 440 

Ctm INITIALIZE DATA FOR TIME INFORMATION BLOCK TT-TT 

AIBT<1> = 0.0 
Du 450 I— 2 .* C 

450 AIBTU :> = AIBK I > 

CTT IF CALIBRATING SET NUMBER OF AUERAGES TO 1 TT 

IF ^ I C AL . Eu . 7o .' fOjFiUE = 1 

cm START SAMPLING DATA POINTS TITTT 

CALL ZERO-: TIME, A I BT> 

DO 4G0 I=1,N0AUE 


j' 



A6l 

A62 


460 

470 


CALL AN I HP •: IFRAME.. IER> 

IF ( ( I FRAME . GE . 0 > . AND . ( I ER . GE . 0 > > GO TO 462 
TYPE 451. .. I ER.. I FRAME 

FORMAT ‘ERROR ANINP ‘ .. 15.. 15>L ' FRAME CODE 

STOP 

IF (ICAL.EQ.78> UJ TO 50O 
HrBD:2.':-=AI&Iv2;:' 

CALL FLO I NF< I ARRAY .. A I B I .. CARR AY .. A I BC 
CALL ADD( CARRAY , A I BC , T I ME , A I BT > 

CONTINUE 

DO 470 N=l.. ISIZE< INDEXZ'.:- 
T I ME( N ;>=T I MEv N >■' NOAUE 


15 > 


CTT CONUERT TO PRESSURE UNITS OF N/(MTM:> TT 

480 BETA=10.0T:p:-OIAIN>/20. .> 

PRESUR= •: BETAr< 31 . Z/ '.'CALv INIK > > > 
CALL MLCuNRv PRESUR , T I ME • A I BT 
Lu= 1 


cmi GRAPH TIME HISTORY mmt 
GO TO 530 

C»;TTTTTTT:1 TTTT l-TTT I- TTTTT 

CT DATA FOR POIJER SPECTRUM GRAPH A 


CTT CONUERT DATA TO NORMALIZED INTEGER FORM T-T 

450 AIBK2>=AIBT(2> 

CALL F I >(F I N( T I ME • A I BT . I ARRAY . A I B I > 



Cl^t: DIRECT FOURIER TRANSFORM TT 

500 CALL DFTv I ARRAV ..AID I .. I ER 1 -■ 

IF vIER:.EQ.0> uu TO 520 
TYPE 510. lER 

5 1 0 FORMAT ( /•••'•• .. 5>^ .. ' ERROR DFT ‘,15 > 

STOP- 

CTT CORRECT FOR ONE-SIDED DFT YT 

520 DO 521 I =2 , 1 S 1 2Ev I NDEXZ > 

52 1 I ARRAVv I >= I ARRAY-: I >T2 


CTT OLD CALIBRATION ? THEN NO NEED TO INITIALIZE DATA FOP 
CTT TIME INFORMATION BLOCK TT 

■ i' 

IF •: IChL.EQ.79> go to 52S 

CTT CALIBRATING SO INITIALIZE DATA FOR TIME INFORMATION BLOCK 

AIBT-:i> = 0.0 
AIBT(2:-=AIBI(2;.- 
DO 525 1=^3, e 

525 hIBT-: I ;-=AIBI-: I > 

526 CALL ZERO •: TIME. AIBT > 

AIBT(5:>=4.0 

m PERFORM SELF-CONJUGATE MULTIPLY-AND-ADD FOR AUERAGED 
CTT AUTO SPECTRUM FROM DFT RESULTS YTT 


CALL ASPEC-: I ARRAY .. A I B I , T I ME , A I BT > 



CONUERT DATA TO FORTRAN FLOATING POINT FORM ^4: 

CALL F LOTFC T I ME .. A I BT > 

C44c:r CORRECT UALUES TO RMS PRESSURE SQUARED 

CALL MLCUNRC . 5 , T I ME . A I BT > 

CTT IF CALIBRATING THEN SUM MAGNITUDES SQUARED 
CTT FOR CALIBRATED DATA m 

IF ( ICAL.EQ.7y> GO TO 540 

C4:4: IF HOT CALIBRATING THEN GROUP POHER SPECTRA IN 2‘S tt 

I MAX. = ISIZEaHDEXZ>/2 
DO 520 I = l.IMAX 
IPI = I+I 
IPIMI = IPI - I 

528 TIME-; ! • = TIME< IPI > + TIMEdPIMI ) 

CTT SET LONER LIMIT FOR PuNER SPECTRA TT 
DO 530 1 = 1.. ISI2EnNDE>:Z> 

530 IF •: T I ME( I > . LT . < 4 . OF-'. 1 0 . i:+:( -y . > > > > T I ME( I >=4 . mx 1 0 . tU 
m. GRAPH POIIER SPECTRUM m 


LG = 2 
GO TO 630 

Ctm SUM MAGNITUDES SQUARED FOR CALIBRATED DATA tm' 


540 


UCAL< IMIK>=0 



DO 55ti 1=400.600 

550 UChLv I H I K >=UChL< 1 1 1 1 K >i-T 1 ME*- I > 

UCALv IMIK >=SQF:;T< UChU IHIK > > 

TYPE 560,. I NIK 

560 FORNOK ‘WUHr IS THE ChLIBRATIuM FhCTuR FDR NICROPHOME 
1 13..* vREAL IMPUT :>‘ > 

ACCEPT 570, CAL.FhC 
570 FORMAT v 1 8 . 3 > 

UCAL< IMIK >=UCAL< IMIK>TCALFAC 
TYPE 580 . 1 M I K . UC AL< I M I K > 

5Sy FuRMAT(8K/ ' UChLv ‘ 12.. ‘ > = * ^F8 6 > 

CKIKK PROMPT FOR NEXT MICROPHONE TrTITTT 


GO TO 110 


CTTTTTTTTT r TTT TKK K 

LT DATA FOR SOUND PRESSURE LEUEL GRAPH I 
TKKTKKKTTTTK TIT 


CTT CuNUERT TO DB SCALE TTTT 

600 A = 2.513 10 . TT0. > 

C ALL MLCONR-: A , T I ME , A I BT > 

DO 610 I = I .. I S I ZEv I NDEX2 > 

610 IF < TIHEv I ;> LT lO . > TIME< I >=10 
CALL LUG(TIME,AIBT> 

C ALL MLCOHRi I O . , T I ME , A I BT • 

00 602 I = 1 , I S 1 2£< I NDEX2 >' 2 
TIME<i ;> = TIME-: I :>-50. 

602 IF < TIME( I > . LT . O > TIMEX I >= O . 

LG=3 



C^: GRAPHS t 


Ctm DRAW BOhiRDERS rmi- 




G40 


650 


CALL ERASED ‘ KB ‘ > 

CALL BEAMP< ' KB ‘ . 1.41. .• 122 > 

CALL UECTuRC 'KB' . 102S.. 122:> 

CALL UECTORC 'KB' 1023..?20:> 

CALL UECTuRv ' KB '.141.. 720 > 

CALL UEC;T0R< ' K!B ' .. 1 4 1 .. 1 22 > 

IF V LG EQ . 1 :> CALL BEAMP':; ' KB ',.141.421 > 

IF ( LG . ED . 1 > CALL UECTORv ' KB' ,. 10i^3. 421 > 

IF ( LG . EQ . 1 > CALL ALPHA-'. ' KB ' . 460. 760,. 12.. 

IF LG , EQ . 2 CALL ALPHA-' ' KB ', .,442. 760,. 14. 

IF -:LG EQ.?.:- CALL ALPHA-', ‘ KB 392, 760. 20. 
EMCODE-: 1 2. 640. AC HAD > HOAUE 
FORMAT-:. 13. IX. 'AUER AGES' :• 


TIME HISTORY' j- 
POWER SPECTRUM ' 

SOUND PRESSURE LEUEL' 


CALL ALPHA-.: ' KB ' . 4 60. 729, 1 2. ACHAU > 


IF •:lg.eq.i> 

CAL.L 

ALPHA-: ‘ KB ' 

. 1 . 500, 9. ' 

AMPLITUDE' 

IF (LG.EQ.i;:- 

CALL 

ALPHCi-:. ' l:b ' 

. 1.461.7.. ' 

-: N. rvKM ;> ' ;:- 

IF (LG EQ.?:- 

CALL 

ALPHA-: ' KB ' 

. 35. 407. 2. 

' DB ' :> 

IF (LG.EQ.3:- 

CALL 

AL.PHA-: ' KB ’ 

. 30. 368. 3. 

• REF ' ;- 

IF (LG.EQ.?.:- 

CALL 

ALPHA-: ' KB ■ 

c* 1 

00002 ' ’> 

IF (LG.EQ.2> 

CALL 

ALPHA-: ' KB ' 

> 1 • o/ * 

PRESSURE ' ;:- 

IF (LG EQ.2:- 

CAL.L 

ALPHA-: ’ KB ' 

. 1.461.7. ' 

SQUARED' > 

IF (LG.EQ.2;- 

CAL.L 

ALPITA-: ‘ i:b ' 

. y 7. 422. 1 .. 

' 2 ' > 

IF ^ LG. EQ.2 > 

CAL.L 

ALPHA-: 'KB' 

. 1. . 399. 7. ' 

-: Tt - mn :> ' ;:- 

IF (LG.GT.l :- 

CAL.L 

ALPHA-: ' KB ' 

.471. 50 , 1 4 

. ' FREQUENCY •: HZ ;> 

IF (LG .EQ 1> 

CALL 

ALPHAv ' KB • 

. 500. 50. 10 

. 'TIME -:SEC;-' 

ENCODE-: 60. 65( 

D . AC HAR ;:- I TEST . I RUN . I PO I NT 

. IMIK 

FORMAT-: * TEST 

' . I 

3. 4>: . ' RUN 

il 

'POINT = '.13.. 


1 4X .'MI CROF HONE = '.12. 3X 

CALL ALPHA( ' KB ' . 1 45. 20. 60. ACHAR > 



Y-hXIS scale 


C'y 


IF <LC HE.l> H = ISIZEv INDEKZ V2 

IF (Li; EQ.r> N = isize(inde:>:z> 

I F *' •' LG . Em . 1 > . AMD . < M . GE . 5DD > > N = 500 

CrU. FIMD NhXIMUM AMPLITUDE mm 

hMAXAP = ABSiTINEa ;:■> 

Du 660 1=2, N 

66u I F ( AMhXhP . LT . ABSi T I HE( I > > > AM AXAP = FiBSC 

IF vLG.EQ.l) GO TO 6S0 


TIME< I >> 


CTTTtT FOR SPECTRA PLOTS DOM'T GRAPH ANYTHING LESS THEM 50 DB tt 
DO 670 I=1,H 

670 IF‘! ! I ME'.. I X LI . 00O04 ) TIME( I ) = .00004 


cmmm. compute scaling factor mmmt 

68£i TEST=ALOG 1 0( AMAXAF' ’> 

TIF-hIMT< test .:- 
FF -TEST-TIP 
FACT = 2. 

IF <FP GE. 25:;- FACT = 5 
IF ( FP IjE 66 ) FACT = 10 
IF (FP.GE. .04> FACT = 20 . 
SCALEY=< 1 0 . m.T I P ;-TFACT 
r3CL.Y=INT-; SCALEY > 

IF (LG.EGLS) GO TO 7'20 



t 


TIC MARIAS FOR V AXIS GRAPHS I AND 2 TTTTTT 

ITICV = 122 
DO 6o5 1=1.10 

CALL- BEANP-:; * m * • H 1 . I T ICY > 

CALL UECTORv ‘ KB M 25 . H I C Y > 
c-85 ITICY = IT ICY + C0 

CALL BEAMPv ‘ l:B * 1 4 1 . 72 1 • 

CALL DECTORv ‘KB* 125..?21 > 

CTTTTTTT LABELS FOR Y AXIS TIC MARKS FOR GRAPHS 1 AND 2 TTTT-T 

IF ( LG . EC! . 2 > CALL ALFHA( ' KB ' , 60.. 122 / 2 ^ ' O . ' > 

IF -'LG ME . 2 > ENCODE-: 5 .. GSO ,. ACMAX > I SCL.Y 
6S»0 FORMAT •: 14.. ‘ ‘ > 

IF ■: Ll . NE 2 ■' CAL-L AL PHAC * l:B ' ' > 2o .. 72o .■ Li ACMAX > 

IF LL . EQ 2 ' ENCuDE*: ‘2.. 700 . AL2MAX ISL^L. r* 

700 FORMAT':; IS.. ‘ ‘ ;• 

IF •: U.j . EO . 2 ;• CALL ALPHA': ' KB ' .. 1 .. 720 . 9 - AC2MAX ;• 

IF •: Ll EQ . 1 ;• CALL ALPHA*: ‘ KB ‘ / 60.. 421.2 . ‘ 0 . ' 

IF •: LG ECi . 1 ENCODE': 6 .. 7 1 C> .. ACM I N :• I SCLY 
710 FORMAT': ‘-V. 14. 

IF ( LG^EQ . 1 CAL.L ALPHA*: ' KB ' / 20.. 122.. 6 . ACM IN .> 

GO TO 740 

mm. TIC MARfxS FOR Y AXIS GRAPH 3 mm 

720 ITICY=122 

DO 730 1=1 /S 

CALL BEAMP< ‘ KB ‘ , 1 4 1 .. I T I C Y > 

CALL UECTORC ‘KB* .. 125. ITICY ;• 

730 ITICY=ITICY+75 

CALL BEAMP*: * KB ‘..141.. 72 1 

CALL UECTOR': ‘KB* / 125. 721 > 



LABELS FDR V A>as TIC MARKS GRAPH 3 ^ 
CALL ALPHA-; ‘ KB ‘ 4D, 72D. 4 . ' 130 . ’ 



CALL ALPHAC ' 

KB‘ ,.55.. 

122. .3. • 50 . ■ 


cm. COMPUTE X AND V 

AXIS 3 

CALING FACTOR FOR 


740 IF (LG.EQ.1> 

SCTiLEX 

= 1 . 


IF •;';lg.eq.i 

> FiND . •; 

I size:-; indexz '.llt. 


IF <( LG .EQ. 1 

> . AND . ( 

isiZE'; indexz:>.eq. 


IF vLG.HE. 1 .> 

AMAXAP 

= SAM-; lNDEXS>/2 


IF (LG.NE.l) 

3CALE>-! 

= 17d4./IMAX 


HP = INT*-. oc'2 

'SCTiLEX 

'i 

•' 1 . 


IF (LG.EQ l> 

VOFSET 

= 421 ' ' 


IF vLG.EQ.l) 

SCALEV 

= 2^9 . ^ SCALEV 


IF vLG.GT.i;:' 

VOFSET 

= 122. 


IF (LG.EQ.2> 

SC 1 iLEV 

= 5S'S . ^'SCALEV 


IF <LG,EQ.3> 

SCCiLEV 

= 52S . /'SO . 


SCALE>; 

SChLE>=: 


cm. TIC HARKS FOR >=: Ar-ilS TTTTTT4;TT4: 

IF'.; LG . Eu . 1 > hHAXAP=NP/( eDHTISIZE-; INDEX?. > > 

test=aluli0'; ahaxap > 

TIP=AINT( TEST;- 
FP=TEST-TIP 
IF (FP LT.0;> TIP=TIP-1 
IF <TIP.LT.0> FF-FP+1 
FACT=I . 

IF < FP . GT . 0 . 30 1 02 > > FACT = 2 . 

IF ^ FP . G1 . 0 . CAADF ;• ’> FACT 5. 


1 



T Ac. 


?44 

750 


UhLKEV == FACT T vlA.TTdlP):' 

0ELNRK=UALKEV^' 10 . 

MMARKS=1 + INK hNhXhP' DELMRK > 

DO 750 H= 1 .. NMhRKS 
MML^N-l 

I XMRK = 141 + 1 NT ( NnLTDELriRFT0:32--HHHXHp > 

CALL BEhMP-: ‘KB' • IXNRK.. 122.:' 

CALL OECTORv ‘KB‘ LCMRK. 1 lO;' 

IF •: N , EC! . 1 CALL ALPHA-: ‘ KB ' .. IKMRK.. 75, 2.. ’ O ‘ > 

IF 'rN HE ll;* GO TO 750 
CALL BEAMP( ‘K.B‘ .. IKMRK,. 122 > 

CALL UECTOR': ‘KB' , IKMRK, 1O0> 

IF *: LG . EO . 1 ENCUDE 7.- 742.- AcTMAR OALKEV 
F0RMHr':F7.5:.' 

IF •: LG . EQ . 1 > CALL ALPHA'.' ‘ KB ’ .. IXMRK-TO.- 75.. 7.- ACTMAR> 
IF LG . HE . 1 > ENCODE': £• • 7' 44.. ACHAX > '.'ALKEV 
FuRHhT': FC7 0 > 

IF ■: LG , NE . 1 ' CALL ALPHA-: ' KB ‘ .. IXMRK-TO, 75, C.. ACNAK > 
CONTINUE 


CTTT PLOT POINTS. FOR ALL 3 GRAPHS 

IV = I NT <v UPSET + TIME'.: 1 ;:':TSC ALEV 
IX = 141 . 

CALL BEAMP-: ‘KB' .. IX, IV ;• 

DO 760 1=2, HP 

IX= I NT-: 141. + -: I-1>:TSCALEX .:' 

I V= I NT'.: VOFSET+T I HE-: I .>TSCALEV > 

760 CALL UECTOR-: ' KB * , I X , I V > 

CTTTTT DONE HITH PLOTTING NOW NAIT FOR COMMAND TO CONTINUE 

ACCEPT 770 
770 F0RMAT<1A1> 



Ctm IF DONE WITH TIME HISTORY THEN GO hND ChLCULATE DOTh FOR 
crm POWER SPECTRUM GRAPH ll-TT 


IF (LG.EQ.1> GO TO 45*0 


CTTT IF DONE WITH POWER SPECTRUM GRAPH THEN GO AND CALCULATE 
cm DATA FOR SOUND PRESSURE LEUEL GRAPH TTTT 

IF (LG,EO,2;> GO TO GO0 

CTT IF DONE WITH SOUND PRESSURE LEUEL GRAPH CLEAR 
CTTT SCREEN AND START OL'ER m 

CALL ERASEv * KB ‘ > 

GO TO 323 

CT PROcRAM ENDS X 

S9S STOP 
END 

cxxxxxxx 

c 

C ECHO SUBROUTINE, TO ECHO CHARACTERS ON SCREEN 
C 

SUBROUTINE ECHO •:iCAL> 

CALL I POKE < "44, "00000 > 

I=ITTOUR( ICAL> 

CALL I POKE •: "44, "10000 ;> 

RETURN 

END 



APPENDIX D 


EXAMPLES OF TYPICAL RUNNING SEQUENCES 


D-1 



ENTER TEST NUMBER NNN 1 
ENTER RUr^ NUMBER NNN 1 


ENTER POINT NUMBER NNN 1 

NHhT type of CHL.IBP.HTIOri ? 0 OR N. vNO <CR>> N 
ENTER NUMBER OF MICROPNONE TO ChLIBRhTE 01 

uhht is the ChlibrhTion Fhctor for microphone 
1 . 

UCHL*' 1 > = 0 . 57yyS4 

ENTER NilMBER OF MICROFl-lONE TO ChLIBRhTE OS 

HHhT is THE ChLIBRhTIOM FhCTOR FOR MICROPHONE 

I . 

UPhL'- > = 0 5?yOr.l 

ENTER NUMBER OF MICROPhlONE TO ChLIBRhTE 07 
MHOT IS THE ChLIBRhTIOM FhCTOR FOR MICROPHONE 

1 

[)' hL'- t" ‘ ~ 0 . SyuOO.!' 

ENTEFr NUMBER 0^ MICROPHONE TO ChLIBROTE 00 


1 <REhL INPUT 


5 CREhL input 


7 <REhL input 



I'OLTS e:hte:r shnF'LE--freg!Ue:iicv cut-off frequeucy enter 




u:> r-- 


Tvj r-vj r^j r-^j 




i^} 

cvj »::5 

— »>j 




^ AJ 


UJ 

o 


LU 

O 

Ci 

i‘ ’t 









u 

1 > • 


01' 

'Hi 

‘H 

'H 



LlJ 

Q' 




. 





UJ 





'H* 



''’X 

Ct: 

b" 



T7 



‘H 

h- 

Ll. 


A] 

u:i 




*‘=J 

! 

1 




'Aj 

U^' 

'H 



LU 






AJ 

U’ 


-J 









!i 








Cii 









LU 

<X 









Lf' 








Lu 

LL 









LU 


















”7“ 









LU 



r^:j 

tT 

U" 

^JL» 





m 

'IT* O 

^ OJ rt OD 
) i } 1 1 } i 


+ -f + 4- + + 


\ 


0-3 



HUMBER OF ShMRI.es 

1 

5 IS 
1024 
.^0*1 y 

ENTER SAMPLE CuDE/2 
BANDWIDTH = lij.uyA 
UOL.TS = 4 L10 

SAMPLE FREu.UENCV = 5120. 0U 

ENTER MICROPHONE I IN 0? 


ENTER 


ENTER CAIN (REAL INPUT > O. 

GAIN IS 0.0 V OR N (NO <CR>;:. i 

DATA ACQUISTION STARTS WHEN <CR> ] 
ENTER NUMBER OF Al'ERAGES NNN 00 1 


S ENTERED 






OUIICI F'RESSURE L.EUEL 
I h'.,.'ERhGES 






i£i 

a 


Li- 

LiJ 

0 ^ 


izii 


o 


D-7 


0-3 


ENTER 


NUMBER OF SAMPLES 

1 r.!o 1 

51 ?: 2 

L02H 3 

Su-IS 4 

ENTER SAMPLE CODE ,i2 

BhNDN I DTH = 1 Cj . 000 

UOL.TS = +• - 4.00 

SAMPLE FREOUENCV = 512O.00 

ENTER MICROPMONE NN 07 

ENTER GAIN (REAL INPUT > O. 

GAIN IS 0.0 V OR N (NO <CR>> 

DATA ACQUrSTIuN STARTS NHEN <CR> 
ENTER NUMBER OF AUERAGES NNN Olo 


IS ENTERED 






MICRuPHuNE 




10 hUERhGES 



POINT = 1 MICROPHONE 





HUMBER OF SAMPLES EHTER 


O 

UJ 

a: 

111 

UJ 

CO 


Ci:' 

ill »U) 





0-12 


EHTER HLlMeER OF AUERAGES HHH OlO 



1 0 hUE:Rh*.^E 
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